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ABSTRACT 

Considered  herein  are  the  partial  differential  equations, 

u^  +  Ujj  +  uux  +  Lu  *  0,  (*) 

3  2 

where  L  denotes  3^  (equation  A)  or  -3^3^  (equation  B).  In  (*)  u  is  a 

real- valued  function  defined  for  all  real  x  and  for  t  >  0,  and  interest 
will  be  focused  on  solutions  of  (A)  and  (B)  that  correspond  to  the  initial 
condition  that  u(x,0)  is  a  given  function*  Equation  (A)  is  the  Korteweg-de 
Vries  equation  (1895)  while  (B)  is  the  model  studied,  for  example,  in 
Benjamin,  Bona  and  Mahony  (1972).  It  has  been  argued  that  either  (A)  or  (B) 
can  be  used  with  equal  justification  to  model  various  physical  phenomena. 

To  establish  this  claim  an  exact  relation  connecting  solutions  of  (A)  and 
(B)  is  derived,  showing  that  the  two  models  yield  predictions  whose  difference, 
over  significant  time  scales,  is  only  of  the  small  order  that  is  formally 
neglected  by  either  model. 

Complementing  the  theoretical  study  are  some  numerical  experiments  based  on 
(B).  These  experiments  suggest  that  the  aforementioned  theoretical  estimates 
are  sharp,  and  that  they  are  valid  *v  to  the  time  scale  for  which  either 
equation  formally  ceases  to  be  an  accurate  model  of  underlying  physical  phe¬ 
nomena.  The  experiments  also  indicate  that  (B)  has  the  property,  which  is 
well  known  for  (A),  that  certain  classes  of  initial  data  evolve  into  a 
sequence  of  solitary  waves  followed  by  a  dispersive  wave  train. 
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1.  INTRODUCTION  1St  Spcc,ial 


-"This  paper  is  concerned  with  mathematical  models  representing  the  uni¬ 
directional  propagation  of  weakly  nonlinear  dispersive  waves.  Interest  will 
be  directed  toward  two  particular  models  that  were  originally  studied  in  the 
context  of  surface-wave  phenomena  in  open-channel  flows. 1  The  rationale  behind 
the  derivation  of  these  models,  details  of  their  mathematical  properties  and 
their  applicability  to  a  host  of  quite  disparate  physical  systems  are  well 
documented  (see,  for  example,  the  review  articles  of  Benjamin  1974,  Bona  1980, 
1981,  Jeffrey  and  Kakutani  1972,  Kruskal  1975,  Miura  1974,  1976,  Scott,  Chu 
and  Kf'Laughlln  1973  and  Zabusky  1981). 

In  attempting  to  describe  open-channel  flows  the  underlying  principle  in 
the  derivation  of  these  model  equations  is  that  their  solutions  should  approx¬ 
imate  solutions  of  the  two-dimensional  Euler  equations,  posed  with  appropriate 
boundary  conditions  at  the  bed  of  the  channel  and  at  the  free  surface.  Within 
the  context  of  procedures  for  generating  such  models  it  is  possible  that 
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several  different  equations  may  emerge.  The  choice  of  which  approximation  to 
use  will  then  depend  on  properties  of  one  model  vis-a-vis  those  of  another. 
For  long  waves  on  the  surface  of  water,  two  models  have  received  particular 
attention.  One  is  the  equation  of  Korteweg  and  de  Vries  (1895)  (equation  (A) 
or  the  KdY  equation  hereafter)  and  the  other  is  an  equation  first  studied 
theoretically  by  Benjamin,  Bona  and  Mahony  (1972)  (equation  (B)  hereafter). 

The  qualitative  mathematical  properties  of  solutions  of  these  two  models  have 
been  studied  in  detail.  This  theory  is  rich  and  interesting,  but  appears  to 
offer  no  definitive  reason  for  preferring  one  or  the  other  of  these  models 
for  the  task  for  which  they  were  originally  derived. 

.''The  purpose  of  the  present  paper  is  to  make  a  quantitative  comparison 
between  the  solutions  to  tne  initial -value  problem  for  each  of  these  models. 
The  basic  conclusion  of  the  study  is  that,  on  a  long  time  scale  T  naturally 
related  to  the  underlying  physical  situation,  the  equations  predict  the  same 
outcome  to  within  their  implied  order  of  accuracy.  In  this  case  the  choice  of 
one  of  these  nodels  over  the  other  to  describe  a  physical  problem  is 
apparently  immaterial,  with  factors  of  incidental  convenience  probably  provid¬ 
ing  the  main  criteria  in  a  given  situation.^  It  is  worth  noting  that  if  one  is 
only  interested  in  the  evolution  over  a  much  shorter  time  interval  than  [0,T], 
then  a  model  simpler  than  either  of  the  aforementioned  (a  factored  version  of 
the  linear  wave  equation)  will  suffice.  The  main  analysis  leading  to  the 
above-stated  conclusions  is  presented  in  section  four.  The  earlier  sections 
two  and  three  give,  respectively,  a  brief  account  of  the  assumptions  and 
formal  limitations  Inherent  in  the  models,  and  some  mathematical  definitions 
and  results  needed  for  the  analysis  of  section  four. 

The  question  of  the  relationship  between  the  two  model  equations  has  been 
discussed  in  general  terms  by  Benjamin  et  al.  (1972),  Kruskal  (1975)  and 
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Whltham  (1974).  Showalter  (1977),  working  carefully  through  the  standard 
formal  scalings  and  expansions  leading  to  models  such  as  those  considered 
here,  derived  some  alternative  systems  and  conjectured  explicitly  that  the  KdV 
equation  and  equation  (B)  will  give  similar  answers  on  a  time  scale  much 
smaller  than  the  scale  T  mentioned  earlier.  Bona  and  Smith  (1975),  in  their 
paper  on  the  initial -value  problem  for  the  KdV  equation,  deal  with  exactly  the 
issue  considered  herein,  but  give  no  attention  to  the  time-scales  over  which 
their  results  are  valid.  For  practical  purposes,  these  time  scales  are 
crucially  Important. 

In  relating  the  solutions  of  the  two  model  equations,  the  present  work 
supplements  the  studies  of  Berger  (1974)  and  Kano  and  Nishida  (1953)  in  which 
the  relation  of  the  model  equations  to  the  Euler  equation  is  examined.  It  is 
also  complementary  to  the  various  studies  comparing  predictions  of  these  model 
equations  with  the  outcome  of  some  laboratory  experiments  (see,  for  example. 
Bona,  Pritchard  and  Scott  1981,  Hammack  1973,  Hammack  and  Segur  1974,  and 
Zabusky  and  Galvin  1971). 

In  addition  to  the  theoretical  relation  linking  the  two  equations,  we 
have  also  made  numerical  experiments  designed  to  afford  further  comparison 
between  the  two  models.  The  first  set  of  numerical  experiments  relates 
directly  to  the  results  of  section  four,  showing  how  particular  solutions  of 
the  two  models  differ  as  a  function  of  time.  The  second  experiment  appears  to 
confirm  that  equation  (B)  shares  the  property  with  the  KdV  equation  whereby 
certain  classes  of  initial  data  evolve  into  a  sequence  of  solitary  waves 
followed  by  a  dispersive  tail.  These  results  are  reported  in  section  five. 


2.  THE  MODEL  EQUATIONS 


A  model  often  used  to  describe  the  unidirectional  propagation  of  irrota- 
tional,  weakly  nonlinear, dispersive  waves  on  the  surface  of  an  ideal  liquid  in 
a  uniform  channel  Is  the  equation  proposed  by  Korteweg  and  de  Vries  (1895), 

"t  *  \  *  7  n\  +  £  "xxx  *  °‘  ,,a) 

In  this  equation  n  *  n(x,t)  represents  the  vertical  displacement  of  the 
surface  of  the  liquid  from  its  equilibrium  position,  t  is  the  time  and  x  Is 
the  horizontal  coordinate  (which  Increases  in  the  direction  of  propagation  of 
the  waves).  Equation  (la)  Is  written  In  dimensionless  form,  with  the  length 
scale  taken  to  be  the  undisturbed  depth  h  of  the  liquid  and  the  time  scale  to 
be  (h/g)1/2;  g  is  the  gravity  constant. 

It  Is  assumed  in  the  derivation  of  (la)  that  the  maximum  amplitude  e  of 
the  waves  Is  small  and  that  the  waves  can  be  characterized  by  a  horizontal 
scale  6"1,  which  Is  large.  In  particular.  It  is  crucial  that  the  amplitude 
scale  and  the  horizontal  scale  of  the  waves  are  such  that  efi”2  Is  of  order  one 
so  that  the  nonlinear  and  dispersive  corrections  to  the  primary  terms  *nd 
nx  are  of  comparable  importance  (see,  for  example,  Meyer  1979).  These 
considerations  suggest  the  Introduction  of  a  new  dependent  variable  N  and  new 
Independent  variables  $  and  t  such  that 

n  3  cN,  x  3  e  t  3  e  (2) 

in  which  case  N  and  Its  derivatives  are,  by  assumption,  all  order  one  quantities. 
In  these  variables,  equation  (la)  may  be  rewritten  as 

Nt  +  NE  +  7  +  F  eNc«  “  0(e2),  (lb) 

where  the  order  of  the  terms  neglected  in  the  derivation  of  (1)  has  been 
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Indlcated  explicitly  on  the  right-hand  side.  It  is  apparent  from  (lb)  that 

3  1 

the  nonlinear  and  dispersive  terms,  j  and^-eN^,  resPect*ve1y»  con" 
stitute  corrections  of  order  e  to  the  order  one  primary  terns.  Thus  a  formal 
calculation  based  on  the  simple  equation. 


suggests  that,  on  the  time  scale  *  e"1,  the  nonlinear  and  dispersive  terms 
may  have  had  a  significant  Influence  on  the  structure  of  the  waves.  In  the 
coordinates  appearing  In  (la),  corresponds  to  a  time  scale  t}  =  e"3^2.  By 
the  same  reasoning  the  terms  neglected  In  the  derivation  of  (1)  could  have  had 
a  cumulative  effect  of  order  one  on  the  time  scale  t2  s  e'2  (corresponding  to 
the  time  scale  t2  *  e"®^2).  These  arguments  further  suggest  that  the  equation 
Nt  +  *  0  wuld  suffice  to  describe  wave  evolution  on  a  time  scale  tq  s  1 

(or  t0  *  e-1/2). 

Because  of  the  relative  sizes  of  the  terms  In  (1),  It  has  been  argued 
(see,  for  example.  Peregrine  1966  and  Benjamin  et_  al_.  1972)  that  the  equation 

Nt  +  NC  +  7  cHHz  "  ‘6’  eN«x  “  °»  (3a) 

or  equivalently, 

"t  *  "i  *  7  '  F  "xxt  "  °*  l3b) 

provides  a  model  comparable  with  (1)  for  the  physical  problem  In  question. 

(It  Is  worth  remarking  that  the  present  discussion  applies  to  situations  other 
than  surface  waves  on  an  Ideal  fluid.  Indeed,  one  or  the  other  of  these  equa¬ 
tions  has  been  derived  as  a  model  In  a  wide  range  of  physical  contexts. 
Discussions  of  certain  applications  may  be  found,  for  example,  in  the 
review  articles  of  Jeffrey  and  Kakutanl  1972,  and  Scott,  Chu  and  McLaughl1n 


1973  as  well  as  in  the  text  of  Lamb  1980;  an  account  of  the  principles  under¬ 
lying  the  frequent  appearance  of  these  model  equations  is  given  in  the  article 
by  Benjamin  1974.) 

In  this  paper  we  shall  concentrate  on  the  initial -value  problems  posed  by 
(1)  and  (3)  with 

n(x,0)  =  n0(x),  (4) 

for  xeR*  the  real  numbers.  This  corresponds  to  the  presumption  that  the  wave 
profile  is  known  everywhere  at  some  given  Instant  of  time,  and  that  inquiry  is 
directed  to  the  subsequent  evolution  of  the  wave  field.  Of  particular 
interest  will  be  the  comparison  of  solutions  of  the  two  model  equations, 
subject  to  the  same  specification  (4),  over  the  time  scales  t0,  tj  and  t2* 

To  illustrate  the  kind  of  results  we  have  in  mind,  consider  the  initial - 
value  problems  for  the  linear  versions  of  (1)  and  (3), 


N  +  N  +  i  eN  *  0  (5a) 

t  c  o  see 

and 

•VV  5-  *"«t  ■  °-  (5b> 

together  with  the  initial  condition 

NU.O)  =  MU.O)  =  FU),  (5c) 

where  F  is  order  one.  The  problems  (5a-5c)  and  (5b-5c)  are  easily  solved  by 
taking  the  Fourier  transform  In  the  5  variable.  Let  n,m  and  f,  respectively 
denote  the  Fourier  transform  with  regard  to  5  of  the  functions  M,  N,  and  F, 
the  transformed  Independent  variable  being  denoted  by  k.  We  see  at  once  that, 

n((c,x)  =  exo(-ik[l  -  k2]i)f(k). 


and 


m(k,x)  =  exp(-ik[l/(l  +  J  kfc)3x )f(k). 

Suppose  for  simplicity  that  f  is  smooth  and  has  bounded  support  (f  =  0  outside 
a  bounded  region).  If  the  integral  with  respect  to  k  of  the  absolute  value  of 
the  difference  betweem  m  and  n  is  computed,  as  a  function  of  x  >  0,  there  is 


obtained. 


m(»  ,t)  -  n(*  ,x)  .  -  /  |m(k,x)  -  n(k,x)|dk  <  Ce  t,  (6) 

L1(R)  - 


valid  for  0  <  e  <  1,  say,  where  C  is  a  constant  depending  only  on  f.  It 
follows  that,  for  x  >  0  and  1  >  e  >  0, 


sup  |MU,t)  -  NU,t)|  »  sup  li- /  ei5k[m(k,x)  -  n(k,x)]dk| 


-«<{<« 


-cd<£  <oo 


<  i-  J  lm(k,t)  -  n(k,x) Idk  <  Ce2x, 


where  Z  denotes  C/2*.  Remember  that  M  and  N  are,  like  F,  both  of  order  one 
and  that  the  neglected  order  In  the  variables  In  (2)  is  e.  Then  (7)  shows 
explicitly  at  time  x  of  order  e"*,  when  dispersive  effects  may  have  signifi¬ 
cantly  modified  the  shape  of  the  Initial  wave  profile,  that  M  and  N  are  still 
within  e  of  each  other.  That  Is,  M  and  N  are  the  same  to  the  Inherent 
accuracy  of  either  model  at  time  x  a  e“*. 

Expressing  (7)  In  terms  of  solutions  of  the  linear  versions  of  (la) 
and  (3a)  the  following  result  emerges.  Let  u  and  v  be,  respectively,  solu¬ 
tions  of  the  Initial -value  problems. 


“t  +  “x  *  Suxxx  ■  an<l  \  *  *x  -  F  vxxt  ■  °> 


u(x,0)  «  v(x,0)  =  cF(eA/^x), 
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where  0  <  e  <  1.  Then  there  Is  a  constant  C  independent  of  e  and  t  such  that, 

sup  |u(x,t)  -  v(x,t ) |  <  Ce7/2t,  (9) 

— <x<- 

for  all  t  >  0.  The  inequality  in  (9)  ceases  to  be  interesting  when  t  Is  of 
order  e"^,  since  both  u  and  v  are  of  order  e  in  absolute  magnitude. 

The  principle  object  of  the  present  study  is  to  determine  whether  (9) 
holds  when  nonlinear  effects  are  retained  In  the  model  equations.  More  pre¬ 
cisely,  suppose  g  is  a  given  sufficiently  smooth  function  decaying  appropri¬ 
ately  to  zero  at  t  «».  Let  n  ■  n(x,t;e)  and  c  *  c(x,t;e)  be  the  solutions  of 
the  initial -value  problems. 


3  1 

nt  +  nx  +  7  *x  +  F  nxxx  =  °» 


(10a) 


h  *  5x  *  1  «x  -  k  ‘xxt  X  °* 


(10b) 


n(x,0)  *  ?(x,0)  *  eg(e1/2x). 


(10c) 


CONJECTURE  There  Is  a  constant  C  dependent  only  on  g  such  that  for 
0<e<land0<t<  e"5^2. 


sup  lo(x,t;e)  -  c(x,t;e)l  <  Ztut. 
•<x<» 


In  section  four  this  conjecture,  and  similar  bounds  Involving  derivatives 
of  n  and  c,  will  be  established  for  t  In  the  range  [0,e“3^2].  In  the  next 
section,  notation  and  useful  auxiliary  results  are  set  forth  In  preparation 
for  the  analysis  In  section  four.  In  Section  five,  evidence  from  numerical 
computations  will  be  given  that  the  estimate  is  sharp  and  holds  in  the  full 
range  [0,  e'572]. 
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3.  NOTATION  AND  PRELIMINARY  RESULTS 

The  standard  notation  Lp  =  Lp(R)  will  be  used  for  the  (equivalence 
classes  of)  pth  power  integrate  functions  f:R  ♦  R,  for  1  <  p  <  -,  with  the 
usual  modification  for  p  =  *».  The  norm  of  a  function  f  in  Lp ( R)  is 

if»L  *  {/"  |f(x)|pdx}1/p. 

P  — ’ 

For  non-negative  Integers  k,  H*  Is  the  Sobolev  space  of  L2  functions  whose 
first  k  (generalized)  derivatives  are  also  In  L2.  The  norm  In  H*  Is  taken  to 
be 

.fik  - 

where  f(jj  denotes  the  jth  derivative  of  f.  Of  course,  H°  •  L2  and  we  shall 
use  i  l  to  mean  the  same  as  i  1q.  By  Plancherel's  theorem,  the  norm  in  H* 
may  be  expressed  as  follows, 

ifij;  *  /  (1  +  52k)lfU)l2d£, 

-m 

A 

where  f  again  denotes  the  Fourier  transform  of  f.  If  k  Is  a  negative  Integer, 
then  H*  Is  defined  to  be  the  dual  space  of  H'k.  The  space  Hk  for  k  <  0  may  be 
Identified  with  the  class  of  tempered  distributions  T  whose  Fourier  transform 

A 

T  is  a  Lebesgue  measurable  function  for  which 

»Ti?  a  /  - We)l2d5  <  +  «• 

*  —  1  + 

The  spaces  Hk  for  k  <  0  intervene  only  tangentially  In  our  analysis.  By  If  we 
denote  Ok>0  Hk.  The  elements  of  H"  are  infinitely  differentiable  functions, 
all  of  whose  derivatives  lie  In  L2. 


If  X  is  an  arbitrary  Banach  space,  and  T  >  0,  the  space  C(0,T;X)  is  the 
collection  of  continuous  functions  u:[0,T]  +  X.  If  T  =  +  »,  it  is  required  in 
addition  that  u  be  bounded  for  0  <  t  <  T.  This  collection  is  a  Banach  space 
with  the  norm 

sup  iu(t)iY  , 

0<t<T  * 

where  i  denotes  the  norm  in  X. 

In  the  analysis  given  in  section  four,  the  forms  (la)  and  (3a)  of  the  two 
model  equations  will  be  used.  More  precisely,  by  rescaling  n*  x  and  t  by  the 
order  one  constants  j,  (^■)1^2  and  (^■)1^2,  respectively,  we  may  take  the  model 
equations  In  the  tidy  forms, 

nt  +  nx  +  nnx  +  nxxx  *  (12)  or  (A) 

and 

n^.  +  nx  +  nnx  -  oxx^  *  0.  (13)  or  (B) 

Results  pertaining  to  the  Initial-value  problem  for  both  (A)  and  (B)  will  be 
needed.  Also  intervening  In  our  analysis  is  a  regularized  version  of  (A), 
written  In  moving  coordinates,  namely, 

nt  +  nnx  +  ^xxx  "  enxxt  * 

where  e  >  0.  (In  due  course  the  e  appearing  in  (14)  will  be  identified  with 
the  amplitude  parameter  e  appearing  In  the  last  section.)  The  needed 
theoretical  results  relating  to  (A),  (B),  and  (14)  are  presented  In  the 
following  sequence  of  propositions.  We  commence  with  results  for  the  KdV 
equation  (A). 


PROPOSITION  1.  Let  g  {  Hm  where  m  >  2.  Then  there  exists  a  unique 
function  u  in  C(0,«;Hm)  which  is  a  solution  of  (A)  in  R  x  R+  such  that 
u(*,0)  *  g.  Furthermore,  s£u  C(0,«»;Hm”'^)  for  k  such  that  m  -  3k  >  -1. 

The  correspondence  g  v- is,  for  each  T  >  0,  a  continuous  napping  of  Hm 
into  C(0,T;Hm_3k)  for  ai i  ^  >  q  such  that  m-3k  >  -1. 

Remark.  By  'solution'  we  shall  always  mean  a  solution  u  of  the  differen¬ 
tial  equation  in  the  sense  of  distributions  for  which  the  initial  condition  is 

satisfied  in  the  sense  that  as  t  ♦  0,  u(»  ,t)  -  g(*)  tends  to  zero  in  an  appro¬ 

priate  function  space.  Of  course  If  m  >  3  in  the  above  proposition,  the  solu¬ 
tion  will  In  fact  be  classical.  That  is,  all  the  derivatives  appearing  in  the 
equation  exist  classically  and  are  continuous,  and  the  equation  is  verified 
polntwise  by  u  everywhere  in  the  relevant  domain. 

Proposition  1  summarizes  some  of  the  theory  appearing  in  Kato  (1975, 
1979), and  Bona  and  Smith  (1975).  The  next  result  collects  together  theorems 
2,4  and  5  of  Benjamin  et^_al_.  (1972)  supplemented  by  lemma  2  and  theorem  5  of 
Bona  and  Smith  (1975). 

PROPOSITION  2.  Let  g  Hm  where  m  >  1.  Then  there  exists  a  unique 
function  u  in  C(0,»;H1)  which  is  a  solution  of  (B)  in  R  x  R+  such  that 

u(* ,0)  =  g.  For  each  T  >  0,  u<  C(0,T;Hm)  and,  for  each  k  >  0, 

3^u€C(0,T;H  l).  For  each  T  >  0,  the  correspondence  gi — *u  is  a  con¬ 
tinuous  mapping  of  Hm  to  C(0,T;Hm)  while,  if  k  >  0,  the  correspondence 
g  i — *.3^u  is  a  continuous  mapping  of  Hn  into  CtOJiH"*1). 

In  the  above  proposition,  m  is  an  integer.  Results  of  a  similar  nature 
are  available  for  both  equatir(is  in  Sobolev  spaces  with  non-integer  order 


(see,  for  example,  Kato  (1975,  1979),  Saut  (1975),  Saut  and  Temam  (1976j,  and 
Bona  and  Scott  (1976)  for  equation  (A)  and  Benjamin  and  Bona  (1983)  for 
equation  (B)). 

Both  (A)  and  (B)  have  invariant  functionals  associated  with  the  solutions 
described  in  propositions  1  and  2.  These  will  play  an  important  role  in  the 
subsequent  analysis.  The  simplest  of  these  functionals  corresponds  to  the 
conservation  of  mass.  More  precisely,  if  u  is  a  solution  of  either  (A)  or  (B) 
corresponding  to  initial  data  of  the  kinds  indicated  in  propositions  1  and  2 
(with  the  restriction  for  (A)  that  m  >  3)  then  the  total  ’mass' 

/  u(x,t)dx  (15) 

does  not  change  with  time.  Thus  if  the  Integral  in  (15)  converges  as  an 
Improper  integral  at  t  *  0,  then  It  converges  for  all  subsequent  times  and  the 
value  of  the  integral  is  independent  of  t.  For  (B)  there  are  only  two  further 
invariants  known,  and  Indeed  Olver  (1979)  has  established  that  (B)  has  no 
invariant  functionals  of  the  form 

/  p(u,ux,uxx,...)dx, 

with  p  a  polynomial,  other  than  those  given  In  (15)  and  below  in  (16).  (Such 
Invariants  will  be  referred  to  as  polynomial  invariants.) 

PROPOSITION  3.  Let  g <  Hm  where  m  >  1  and  let  u  be  the  solution  of  (B) 
with  Initial  value  g,  as  guaranteed  by  proposition  2.  Then 

/  Cu2(x,t)  +  ux(x,t)]dx  and  /  [u2(x,t)  -  y  u3(x,t)]dx  (16) 


are  both  Independent  of  t. 


For  the  KdV  equation  (A),  an  infinite  sequence  of  polynomial  invariants 
is  known.  These  take  the  form 

V“>  *  i  [u(k)  '  Cku“(k-U  +  Qk(u’ux . u(k-2>n<,x  ,17) 

where  u^r^  =  3^u  and,  for  each  k  *  0,1,2,...,  is  a  polynomial  of  rank 
k+2.  Here  the  definition  of  rank  employed  by  Miura  et^^l_.  (1968)  is  being 
used.  For  a  monomial,  let 

d  ■  a  p  « 

r<mkiu<o)  “m  —  u(5))  *  j  1 0n  +7>v 

The  rank  of  a  polynomial  Is  then  just  the  maximum  of  the  ranks  of  Its  monomial 
components.  In  fact,  Is  composed  entirely  of  monomials  of  rank  k+2.  The 
next  proposition  is  derived  from  theorem  1  and  proposition  6  of  Bona  and  Smith 
(1975). 

PROPOSITION  4.  Let  gc  Hra  where  m  >  2  and  let  u  be  the  solution  of  (A) 
corresponding  to  the  initial  value  g.  Then  Iq(u),... ,Im(u)  are  Independent  of 
time.  Furthermore,  the  Invariance  of  these  functionals  implies  that  for 
0  <  k  <  m, 

,u(k)(*,t)i  -  qk*,9lk*’  <18) 

independently  of  t  >  0,  where  q^R4  ♦  R+  may  be  taken  to  be  the  square  root  of 
a  polynomial  with  non-negative  coefficients  which  vanishes  at  0. 

Finally,  the  Initial -value  problem  for  the  regularized  equation  (14)  will 
occur  at  a  crucial  point  In  section  four.  For  this  problem  we  have  the 
following  result. 


PROPOSITION  5.  Let  g <  Hm  where  m  >1.  Then  there  is  a  unique  function  u 
in  C(0,»;H*)  which  is  a  solution  of  (14)  in  R  x  R+  such  that  u(*,0)  =  g.  For 

k  m— k 

each  T  >  0  and  integer  k  in  [0,m],  3^u  «C(0,T;H  )  and  the  mapping 

g  | — »  a*^u  is  continuous  from  Hm  to  CCO.T'.H"1"^).  Furthermore,  the  functionals 

V 

j  [u2(x,t)  +  eu2(x,t)]dx  and  /  [u2(x,t)  -^-u3(x,t)]dx  (19) 

are  both  independent  of  t  >  0.  If  m  >  3,  there  is  a  positive  constant  CQt 
dependent  only  on  1913,  such  that  for  0  <  e  <  cQ  and  0  <  t  <  1, 

/  u2x(x,t)dx  <  8q( 1 91 3) •  (20) 

where  aQ:R+  +  R+  is  continuous,  monotone  increasing  and  a0(0)  =  0. 

This  proposition  is  a  minor  modification  of  propositions  2  and  3  in  Bona 
and  Smith  (1975).  Note  that  the  two  Invariant  functionals  in  (19)  are,  apart 
from  the  term  containing  an  e,  the  same  as  for  KdV.  It  follows  as  in 
proposition  4  that  the  solutions  u  provided  in  proposition  5  are  bounded  in 
H1,  independently  of  t  >  0  and  e  >  0. 

COROLLARY  1.  Let  g « Hn  where  m  >  1  and  let  u  be  the  solution  of  (14) 
guaranteed  by  proposition  5.  Then,  for  1  >  e  >  0, 

au(»,t)i  <  igij  and  iux(*,t)i  <  qfig^)  (21) 

for  all  t  >  0,  where  q(0)  =  0  and  q  is  continuous  and  monotone  increasing. 

Proof.  The  first  invariant  in  (19)  shows  that,  for  0  <  e  <1, 
iu(-  ,t)i2  <  /  [u2(x,t)  +  eu2(x,t)]dx  =  /  [g2(x)  +  eg2(x)]dx  <  igi2. 


(22) 


5 


Similarly,  it  is  adduced  from  the  second  invariant  in  (19)  that, 

OD  OD 

IUX(- ,t ) i2  =  igxi2  ♦  J  /  u3(x,t)dx  - -J  /  g3(x)dx 

<  ig«2  +  -j  iu(- ,t)i  ju(.,t)i2  +  -j  ig»2igi^  (23) 

<  ig«2  +  j  igi3  +  j  igi2iu(» ,t)ij, 

where  the  elementary  inequality, 

7  12 

IVI  m  <  I VI i VXI  <  J  I v» x,  (24) 

and  (22)  have  both  been  used.  Adding  (22)  and  (23),  it  follows  that  (21) 
holds,  with,  for  example, 

qU,.[lis2.|s34$4]1/2. 
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4.  ANALYTIC  COMPARISON  OF  THE  TWO  MODEL  EQUATIONS 


Attention  Is  now  focused  on  the  Inltal -value  problems  for  the  model 
equations  (A)  and  (B).  Let  g0  be  a  given,  physically  appropriate.  Initial 
wave  profile.  In  the  scales  Implied  in  both  (A)  and  (B),  g0  must  therefore  be 
of  small  amplitude  and  large  wave  length,  and  these  must  be  appropriately 
related  as  explained  In  section  two.  These  assumptions  on  g0  nay  be  made 
precise  by  Introducing  the  positive  parameter  e,  which  Is  taken  as  a  measure 
of  the  small  amplitude  of  the  wave  profile,  and  assuming  that  gg  may  be 
represented  In  the  form 

gg(x)  «  eg(e1/,2x).  (25) 


here  g  Is  viewed  as  fixed  and  Interest  lies  in  the  regime  e  «  1.  Thus 
consideration  Is  given  to  the  Initial -value  problems 


with 


and 


with 


+  "e"*  *  nxxx  *  °> 

for  (x,t)  cRx  R+, 

i\e(x,0)  •  eg(el^2x). 

for  x  t  R, 

*  'c‘x  -  4t  •  °- 

for  (x,t) e  R  x  R+, 

cE(x.O)  •  cg(c1/2x). 

for  x  c  R. 

(26a) 


(26b) 


In  what  follows,  we  shall  refer  to  a  quantity  of  the  form  q( *  where 
q:R+  ♦  R+  Is  a  C*  function  with  q(0)  =  0,  as  an  order-one  quantity.  Thus  an 
order-one  quantity  Is  a  function  of  a  Sobolev  norm  of  g.  The  principle  result 
of  our  Investigation  may  now  be  stated. 


THEOREM  1.  Let  g  H*  0  where  k  >  0.  Let  e  >  0  and  let  n  and  c  be 
the  unique  solutions,  guaranteed  by  propositions  1  and  2,  of  the  Initial -value 
problems  in  (26).  Then  there  is  an  cq  >  0  and  order-one  quantity  Mj  such 
that  if  0  <  e  <  eg*  then 

inf^f-.t)  -  C(j)(-.t)l  f  MjEj/2  +  7/4 (e 3/2t )  (27) 

at  least  for  0  <  t  <  e"^2,  where  0  <  j  <  k. 

Remark .  We  continue  to  use  the  symbol  U(rj  Introduced  In  (17)  to  denote 
the  rth  derivative  of  u  with  respect  to  the  spatial  variable  x,  a^u. 

4.1  A  Comparison  Over  a  Short  Time  Interval. 

Before  giving  the  proof  of  theorem  1,  a  related  issue  will  be  addressed. 
This  somewhat  simpler  point  is  of  Interest  In  its  own  right,  and  its  resolution 
suggests  an  analysis  of  the  more  complex  situation  reflected  in  theorem  1. 

As  explained  In  section  two,  either  model  equation  (A)  or  (B)  may  be 
viewed  as  a  small  perturbation  of  the  basic  one-way  propagator  ut  +  ux  *  0. 
Moreover,  In  the  formal  derivation  of  these  models  from  more  complete  sets  of 
equations,  assumptions  are  made  concerning  the  sizes  of  various  combinations  of 
the  dependent  variable  and  Its  derivatives.  These  assumptions  play  a  crucial 
role  In  the  derivation  of  the  model  equation,  for  they  allow  certain  terms  In  a 
formal  expansion  to  be  retained  while  others  are  dropped.  As  already  explained 
In  section  two.  It  Is  this  procedure  that  leads  ultimately  to  equations  such  as 
(A)  and  (B).  Moreover,  at  a  cruder  level  of  approximation  than  that  antici¬ 
pated  for  (A)  and  (B),  this  procedure  would  yield  exactly  the  factored  form  of 
the  one-dimensional  linear  wave  equation  ut  ♦  u„  *  0.  It  Is  our  purpose  here 
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to  establish  rigorously  the  natural  suppositions  concerning  the  size  of  the 
dependent  variable  and  its  derivatives.  It  will  also  be  shown  that  either 
model  (A)  or  (B)  may  be  replaced  by  the  equation  u^  +  ux  *  0  without  loss  of 
order  of  accuracy  provided  time  scales  no  longer  than  are  in  question. 

In  view  of  the  theorem  stated  above  concerning  the  relation  between  the 
two  models,  it  suffices  to  consider  only  one  of  these  models  at  the  present 
stage  of  discussion.  Results  established  for  one  model  will  apply  to  the 
other  model  by  virtue  of  theorem  1. 

Consideration  is  therefore  given  to  solutions  of  (26a)  and  their 
relationship  to  solutions  of  the  inital -value  problem, 

o^  +  o^  =  0,  for  (x,t )  iR  x  R+, 

with  (28) 

oe(x,0)  *  eg(e1/2x),  for  x <  R, 


where  0  <  e  «  1.  The  solution  of  (28)  is 


Define 


and 


oe(x,t)  *  eg(e1/2(x  -  t)). 


■ifv  f )  .  — 1  6/  “1/2  .  -3/2*  -3/2* » 

u(x,t)  *  e  n  (e  x  +  e  t,  e  t) 


w(x,t )  =  e“1oe(e'1/2x  +  e'3/2t,  e"3/2t), 


(29) 


(30) 


Then  u  and  w  satisfy  the  Intlal -value  problems, 

u*.  +  uu  +  u  =  ol 
t  x  xxx  1 


for  (x,t)c  R  x  R  , 


(31) 


wt  =  0 


and 


u(x,0)  «  w(x,0)  «  g(x),  for  x<  R. 


Let  h  *  u  -  w.  Then,  h  is  a  solution  of  the  initial-value  problem, 

ht  =  ~uux  '  uxxx*  for  *x,t*  ‘  R  x  R+’ 

and  (32) 

h(x,0)  *  0,  for  x  «  R. 

If  (32)  is  integrated  over  the  temporal  interval  [0,t],  there  appears  the 
f omul  a, 

t 

h(x,t)  *  -  /  [u(x,t )u  (x,t )  +  u  (x,T)]dr. 

Q  A  AAX 

It  then  follows  that  for  any  non-negative  integer  k, 

t 

ih(.,t)ik  <  £  (iu(«,t)ux(.  ,x)ik  +  luxxx(*  .t)ik)dT*  (33) 

This  Inequality  is  the  key  to  the  proof  of  the  following  result. 

THEOREM  2.  Let  g  «Hm  where  m  >  3.  Let  ne  and  oe  be  the  solutions  of 
(26a)  and  (28)  corresponding  to  the  initial  data  eg(c^2x),  where  e  >  0. 

Then  there  are  order-one  constants  Cj  such  that  for  0  <  e  <  1  and  0  <  t, 

,n(d)(" *t)  "  0(j)(,*t)l  -  cjeJ/2  +  9/4t»  (34) 

for  0  <  j  <  m-3. 

Proof.  Fix  e  In  the  range  (0,1]  and  perform  the  change  of  variables 
Indicated  In  (30).  Let  u,w  and  h  be  as  defined  In  (30)  and  just  below  (31). 
For  aqy  T  >  0,  n€ «  C(0,T;Hm).  Hence,  for  any  T  >  0,  u«  C(0,T;Hm).  From  (29) 
It  Is  clear  that  oe  « C(0,T;Hm),  for  any  T  >  0.  Hence  h  C(0,T;Hm),  for  any 


T  >  0.  Thus,  provided  k  <  m-3,  the  right-hand  side  of  (33)  is  finite,  and  the 
formal  calculations  leading  to  this  Inequality  are  easily  justified. 

But,  (33)  Immediately  implies  the  estimate, 


ih(.,t)ik  <  t{iuuyi  k  ♦  »uy  i  k  }, 

K  ‘  x  C(0,t;HK)  xxx  C(0,t;Hx ) 

for  0  <  t  and  0  <  k  <  m-3.  Elementary  considerations  involving  the  Sobolev 


norms  then  yield. 


ih(.,t)ik  <  M.t{iui2  k  .  +  iui  1,4.1). 

*  ‘  *  C(0,t;HK+1)  C(0,t;HK 


holding  for  all  t  >  0  and  0  <  k  <  m-3,  where  the  constants  depend  only  on 
k.  Proposition  4,  which  applies  equally  to  the  KdV  equation  (A)  and  to  the 
KdV  equation  (31)  written  In  coordinates  moving  with  speed  one,  allows  the 
right-hand  side  of  (35)  to  be  bounded  in  terms  of  the  data  g  as  follows: 


ih(.,t)ik  <  Mkt:{q0(igi )fc  +  \+1(igik+1)£  +  Cq0(igi)c  +  V3(l9,k+3l2]1/Z} 


"  Ckt* 


for  0  <  k  <  m-3,  where  the  q^  are  defined  In  (18).  According  to  the  earlier 
definition,  Ck  is  an  order-one  quantity.  It  follows  that 

,h(k)(* •t)l  -  ckt»  (36* 

for  0  <  k  <  m-3.  It  remains  simply  to  express  nc  and  oc  In  terms  of  u  and 
w,  so  Inverting  the  change  of  variables  (30).  This  gives 

ne(x,t)  *  cu(e1/2(x  -  t),e3/2t),  (37) 


and  similarly  for  oe  in  terms  of  w.  It  then  follows  from  (36)  and  (37)  that. 


■  t  /  c  |  xi  >  p  k/2  +  9/4. 

In(kr  "*  ^(k)  *  *  *  1  - 


where  Ck  is  an  order-one  constant.  This  is  just  what  we  set  out  to  prove. 


COROLLARY  2.  Let  g,ne  and  oe  be  as  in  the  statement  of  theorem  2  and 
suppose  M  >  4.  Then,  there  are  order-one  constants  such  that,  for  all 


t  >  0, 


sup  ln/k \(x,t)  -  o^(x,t)|  <  Bktc(k+5)/2 


provided  0  <  k  <  m-4. 

Proof.  First  note  that  if  f «H*,  then  according  to  (24), 

|f(x)|  <  ifi1/2if'l1/2,  (391 

for  all  x  <  R.  Apply  this  with  f  «  and  then  use  (34)  to  bound 

the  resulting  right-hand  side. 


It  is  instructive  to  consider  the  case  k  *  0  in  (38).  The  result  thet< 


reads. 


sup  Ue(x,t)  -  oc(x,t)|  <  B0e2(e1/2t). 
x  <  R  "  u 


Recall  that  both  nc  and  oe  are  of  order  e  In  maximum  magnitude,  and  that 
e2  Is  the  neglected  order  In  the  scaling  appearing  in  (26a)  and  (28).  The 
Inequality  In  (40)  states  therefore  that  n  and  o  agree  to  the  neglected  order 
at  least  over  the  time  Interval  [0,e"^2].  In  other  words,  if  Interest  lies 
In  the  evolution  of  the  wave  profile  over  a  time  t  <  e"1^2,  then  one  might  as 
well  employ  (28)  (that  is,  translate  the  Intlal  wave  profile  at  speed  one 


without  change  of  shape)  instead  of  either  of  the  more  complicated  models  in 
(26).  Over  longer  time  scales,  this  may  not  be  true.  Indeed,  examples  are 
presented  in  section  five  showing  clearly  that  the  estimates  in  (34)  and  (38) 
cannot  be  Improved  in  general.  So  nonlinear  and  dispersive  effects  are 

increasingly  felt  for  times  in  the  range  t  >  e“^.  When  t  is  of  the  order 

the  difference  between  n  and  a  is  of  order  e,  the  basic  size  of  each 

term  separately.  The  models  in  (26)  and  the  simple  model  (28)  have  thus 

diverged  in  their  predictions  by  this  time.  This  accords  with  the  formalism 
explained  in  section  two. 


arison  Over  an  Intermediate  Time  Interval. 


We  now  show  that  solutions  nc  of  the  initial-value  problem  (26a)  scale 
with  respect  to  e  in  the  way  that  is  expected  and  used  in  formal  studies  of 
this  equation. 

THEOREM  3.  Let  g«  Hra  where  m  >  1  and  let  e  >  0.  Let  ne  be  the  solution 

1/2 

of  (26a)  corresponding  to  the  intltal  data  eg(e  x).  Then  there  are  order- 
one  constants  Dj  such  that  for  all  t  >  0, 


_  e  i  ♦  \,  ,  n  j/2  +  3/4 

in( j)  * *t)i  -  * 


whenever  0  <  j  <  m. 


Proof.  This  Is  a  straightforward  consequence  of  proposition  4.  First 
let  u  be  defined  by  (30)  as  before.  Then  u  satisfies  the  initial-value 


problem  (31).  Proposition  4  thus  implies  that  for  all  t  >  0, 


lu(J)<*,t)l  -  qj(l9,j*  s  Dj* 


(42) 
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for  0  <  j  <  m,  as  in  (18).  The  inequality  (41)  now  follows  upon  using  the 
expression  (30)  for  ne  in  terms  of  u  in  (42). 

COROLLARY  3.  Let  g,m,e,  and  ne  be  as  in  theorem  3.  Then,  there  exist 
order-one  constants  Fj  such  that, for  all  t  >  0, 

sup  ln|jj(x,t)l  <  F^e1  +  ^2  (43) 

x<  R 

for  0  <  j  <  m-1. 

Proof.  This  follows  immediately  from  (39)  and  (41)  if  we  define 
Fj  *  <“jVl>1/2- 

Theorem  3  and  its  corollary  show  that  solutions  of  (A)  corresponding  to 
initial  data  of  the  form  eg(e^2x)  scale  with  respect  to  the  parameter  e 
just  as  the  data  does.  In  particular,  (43)  yields 

lnel  -  0(e),  ln*l  -  0(e3/2),  and  lnjxxl  *  0(e5/2), 
as  e+0,  so  that  (44) 

lnen*l  *  0(e3^2), 

as  e+0.  Solving  equation  (A)  for  ti£,  and  taking  account  of  (44),  it  follows 
that 

hjl  ■  0(e3/2), 

as  e+0.  These  results  emphasize  again  that  the  nonlinear  and  dispersive 
terms  in  the  KdV  equation  represent  small  corrections  to  the  basic  wave 
equation  +  nx  *  0.  We  turn  now  to  the  more  exacting  and  technical  task  of 
proving  theorem  1. 


-24- 


Proof  (of  theorem  1).  Suppose  at  the  outset  that  g *  H“,  that  is,  g  is  a 
C“  function  all  of  whose  derivatives  lie  in  L2.  As  in  (30),  define 


and 


u(x,t )  =  e"1ne(e'1/2x  +  e"3/2t,  c‘3/2t), 

v(x,t)  -  .“V(«"1/2x  ♦  e-3'2t,  c-3/2t). 


(45) 


A  short  calculation  shows  that  u  and  v  satisfy  the  initial-value  problems 


ut  +  uux  +  uxxx  ■  °- 


*t  +  *vx  ♦  ’xxx  -  £\xt  1  °-  f 


u ( x,0)  =  v(x,0)  =  g(x).  J 


(46) 


Let  w  =  v  -  u,  so  that  v  *  w  +  u.  Then  w  is  seen  to  satisfy 


wt  + 


"x  + 


w 


XXX 


-  ew 


xxt 


■  cUxxt  -  (uw,x* 


w(x,0)  =  0. 


(47) 


We  continue  to  use  the  notation  f(jj  to  denote  3^f. 

The  task  to  be  accomplished  now  is  the  estimation  of  ||  for 

j  *  0,1,2,..,  .  Because  g«  H®,  both  oe  and  cc,  and  hence  u  and  v,  are  C 
functions  of  x  and  t  all  of  whose  derivatives  are  in  L2  with  respect  to  the 
spatial  variable.  This  fact  justifies  the  following  computations. 

Multiply  (47)  by  W(2j)  and  Integrate  the  result  over  R  and  over  [0,t]. 
After  a  few  Integrations  by  parts,  and  taking  account  of  the  fact  that 
w(x,0)  =  0,  the  following  identity  is  seen  to  hold,  as  in  Bona  and  Smith 


(1975,  eq.  7.9), 


*  o  2  ^  00  12 

/  {"(jjU.t)  +  ew(j+1)(x.t»dx  =  2£  /  w(j)^ut,(j+2)  _  (uw  +7W  1(j+l)}dx<ir 


for  j  ®  0,1,2,...  .  This  relation  will  be  used  repeatedly. 

First,  for  j  *0,  there  appears,  after  two  more  integrations  by  parts. 


.2  .  .  2 


J  (w  +  ewx)dx  “2  ^  /  e(wuxxt)dxdt  '  fQ  J  (uxw  Jdxdi. 


From  this  the  following  inequality  is  derived: 


iwi2  <  f  (2eiwnuxxti  +  iuxiB>iwi2)dT, 


where,  as  before,  i  i  denotes  the  norm  in  L2  and,  throughout  tMs  proof,, 

I  denotes  the  norm  in  L„.  By  a  variant  of  Gronwall's  lemma,  it  follows 
that 

c  c  t  c 

iwi  5  6  <e  1  -  1)  <  etc2e  1  *  M0et,  (50) 

where  C1  and  C2  are  bounds  for  j  *uxi—  and  iuxxti  respectively  and  t  is 
restricted  to  the  range  [0,1].  Using  the  results  of  proposition  4  in  section 
three,  and  using  the  differential  equation,  the  following  estimate  can  be 


made: 


sup  iu  <  sup  ia*(-u  -  uu  li 
t>0  ^  t>0  *  ***  * 

<  sup  (,u,5).  .  .U.JU,,,.  .  3.v.,v,). 


As  remarked  In  (24), 


ifi  <  ifnf  i 
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hence,  in  the  notation  of  (18),  with  the  q^  referred  to  g  of  course,  C2  may  be 
defined  by 

sup  ,uxxt'  -  q5  +  (q0ql)1/2q3  +  3(qiq2)1/2q2  *  C2*  (51) 


It  is  even  easier  to  estimate  a  value  for  Cj*, 


sup  iuyi  <  sup  (iu  nu  i)1/2  <  (q.q~)^2  =  C..  (52) 

t>0  x  t>0  x  xx  "  1 

Since  both  and  C2  are  order-one  quantities,  so  is  Mq.  This  result  Is 

already  interesting,  as  will  appear  shortly.  Further  bounds  lead  to  a  better 

overall  picture,  and  to  L.  estimates. 

Integrating  (48)  by  parts,  in  the  case  j  ■  1,  the  following  relation  is 

derived. 


■  t  * 

J  (w2  +  ew2x)dx  =  2e/  /  (wyu 

—oo  0  -CO 


x  xxxt 


)dxdr 


f  "  -j  o 

^  /  (wx  +  3uww„  +2uvywwv)dxdx. 


(53) 


XX  XX  X' 


The  integrand  on  the  right-hand  side  of  the  latter  equation  may  be  bounded 
above  by 

2t,uxxxt"V  *  +  ^x'-'V*  +  2»V-,W,,V 

<  (4.ux.^  +  IVX.,).WXI2  ♦  e(2luxxxt.  +  aiV'V.'l'V'V’ 


Now  using  the  equation  satisfied  by  u,  and  the  results  of  proposition  4,  we 
may  derive  the  following  estimate,  valid  for  0  <  t  <  1. 
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SUP  {2,UXXXt‘  +  2(,Uw'>Uvvv«)1/2Mnt} 


0<t<l 


XX  XXX 


i  t2,s;<-“ux  -  uxxx»' +  2(vj|1/!ho« 


0<t<l 


,1/2 


1  Z%  +  2(ty*l)i/Cq4  +  8<qiq2)1/2q3  +  6(q2q3)1/4q2  +  2(q2q3) 


,1/2. 


Since  Mq  is  an  order-one  quantity,  C4  is  an  order-one  quantity.  Also,  using 
proposition  4  as  before  and  relying  on  proposition  5  and  Its  corollary  for  the 
bounding  of  ivxi-t  there  appears 


sup  {4iu  1  +  iv  1}  <  4(q.qJ1/2  +  SUp  (iv  11  v  1  )1/2 
0<t<l  x  x  1  ‘  £kt<l  x  xx 

<  4<qiq2)1/2  +  C(qa0)(||g||  3))1/2  *  c3  . 

H 

This  latter  quantity  Is  order-one  also.  Thus  at  least  over  the  Interval 
0  <  t  <  1,  the  Inequality, 

■V2  ~  £  (c3,wx|2  +  eC4'V)dT» 

Is  Implied.  It  follows  Immediately  that, 

C4  c3t  c3 

«V  <  e  jr-  (e  *  -  1)  <  ctC4e  0  -  ctMj,  (54) 

3 

for  0  <  t  <  1.  Since  C3  and  C4  are  order-one  quantities,  Mj  Is  also  an  order- 
one  quantity. 

For  the  case  of  a  general  j,  the  procedure  for  obtaining  a  bound  on 
,w(j)‘  1s  simi1ar  t0  that  followed  above  In  the  cases  j  «  0  and  j  «  1. 

Suppose  Inductively  that  for  j  <  k,  where  k  >  1,  there  have  been  established 
bounds  of  the  form, 


,WU)'  J c  ('Cy'lt  -  »  5 


(55) 
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for  t 1  [0,1],  where  C2 j _i  and  C2 j  are  both  order-one  quantities.  The  goal  now 
is  to  establish  the  same  type  of  bound  for  j  *  k.  To  this  end,  consider  the 
equation  (48)  in  the  case  j  =  k.  Using  Leibnitz'  rule,  (48)  may  be  written  as 


I  (w(k)  *  ew(kH),d’1  *  2  l  t  EW(k)ut,<k+2)d*dT 

—00  J  —00 

t  •  k+1  , 

"  2  0  L  j«0  aj{2'W^+l-j)w(j)  +  w(k+l-j)u{  j))w(k)dxdr  * 

Here  the  aj  are  the  constants  that  appear  in  Leibnitz'  rule.  Separating  the 
top  order  derivatives  and  estimating  the  rest  directly,  we  have 


/_  (w(k)  +  eW(kH)>d*  5  2e|,  "'(k>"“t,(k.2),dl 
t  • 

<ww(k)w(k+l)  +  U¥'(k)''(k+l))dx<h 

(56) 

*  l  l.  aj |w(k)w(k+l-j)w(j) ldxdx 

t  “  k+1 

+  2  l  °jlw(k)u(j)w(k+l-j)|dxdr- 


The  Induction  hypothesis  (55)  assures  us  that  on  the  time  Interval  [0,1] 

HW(J)H  and  are  bounded  by  order-one  constants  if  0  <  j  <  k-1 

and  0  <  1  <  k-2.  By  Proposition  4,  lu ( j ) 1 •  ‘“(j)1-  and»  us^n9  the 
differential  equation,  iu^  ^+2)1  are  boun<le<1  bY  order-one  constants,  for 
0  <  j  <  k+1.  Independently  of  t  >  0.  Also,  by  proposition  5  and  its 
corollary,  for  t  [0,1],  iwi,  iwxi  and  i«xxi  are  bounded  by  order-one 
quantities.  Because  of  (24),  iw  1  and  iwi  are  similarly  bounded.  The 
second  Integral  on  the  right-hand  side  of  (56)  is  equal  to 

l  i  (wx  +  ux >w(k)dxclT  C  l  ,w(k)|ZdT- 
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where,  provided  0  <  t  <  1,  C  may  be  Inferred  to  be  an  order-one  quantity.  The 
third  and  fourth  Integrals  on  the  right  side  of  (56)  are  both  estimated 
similarly: 

t  -  k  *  2 

l  L  jli  *jl“dc)w(kti-j)w(j)l,,xdt :  cl  'W(»‘  * 

*Ul  ljL  'w(k)lMJ,I,)<h 

-C[  '“(k),Z<l1  +eC^  ,‘,(k),dT’ 

valid  for  0  <  t  <  1  at  least.  The  constants  appearing  In  this  Inequality  are 
order-one.  The  same  estimate  holds  for  the  fourth  term  on  the  right  side  of 
(56).  Hence  In  sum,  for  0  <  t  <  1, 

i  lw(ki +  t“(k*nld* :  c2k-i  l  '“(k)'2*  +  eC2k  ( "W*- 

where  C2^_i  and  C2*  are  order-one  quantities.  The  result  (55)  for  j  *  k  now 
follows  and  the  Inductive  step  Is  completed. 

It  Is  worth  noting  that  the  constants  and  C2)c  depend  only  on 

for  0  <  j  <  k+5.  In  fact,  this  consideration  Is  dominated  by  the  term 
ut,(k+2)»  appearing  on  the  right-hand  side  of  (55)  tdilch,  by  use  of  the 
differential  equation  and  proposition  4,  Is  bounded  In  terms  of  qo»...»q|(+5» 
and  so  In  terms  of  igi ••••»|9(|c+5)1  • 

This  last  remark,  coupled  with  the  continuous -dependence  results  In 
propositions  1  and  2,  allows  the  weakening  of  our  Initial  assumption  that 
g  H".  By  approximating  g  Hk+5  by  a  sequence  {gn}^sl  C  if,  we  may  conclude 
that 

■w(J )■  <  ctMj,  (57) 

for  0  <  t  <  1  and  0  <  j  <  k,  where  w  »  u  -v  and  u  and  v  satisfy  (46)  with 
Initial  data  g. 
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Now  it  is  only  necessary  to  translate  the  result  (57),  a  relation 
concerning  u  and  v,  into  a  result  relating  ne  and  ;e.  This  simply  involves 
inverting  the  transformation  involved  in  (45).  It  appears  immediately  that, 

nG(x,t)  =  eu(c1/2(x  -  t),e3/2t),  (58) 

and  similarly  for  and  v. 

Thus,  supposing  now  that  g  Hk+5  and  that  j  <  k,  we  deduce  that 


■  f  Uil1  +  1Cu(j)(£1/2(x-t),e3/2t)  -  ,(j)<£1/2(x-t),e3/2t)])2< 
*  e3+2  ^  <u(j)<z.e3/2t>  -  .(j)(2.e3/2t))2E-1/2dz 

.£3+3/2.wu)(.,e3/2t),2<£3‘7/2M2(c3/2t)2. 

and  this  Is  valid  as  long  as  0  <  e3/,2t  <  1.  Hence  if  0  <  t  <  e"3<^2, 

,fl(J)(*»t)  "  5  ej/Z  +  ?/4Mj(e3/2t).  (59) 

where  Mj  is  an  order-one  quantity.  This  finishes  the  proof  of  the  theorem. 


COROLLARY  4.  Assume  the  hypotheses  and  notation  of  theorem  1.  Then 
there  Is  an  e0  >  0  and  order-one  constants  Nj  such  that  for  0  <  e  <  e0  and  for 
0  <  t  <  e"3/2. 


sup  lnGj)(x,t)  -  C(j)(x,t)l  <  Nje2  +  j/2(e3/2t). 


for  0<  j  <  k-1. 
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Proof.  This  follows  instantly  from  (27)  and  the  inequality  (24)  if  Nj  is 
defined  to  be 

COROLLARY  5.  Let  g«  Hk+5,  where  k  >  0.  Let  e  >  0  and  let  ;e  be  the 
solution  of  (26b)  corresponding  to  the  initial  data  eg(c^2x).  Let  oe  be  the 
solution  of  (28)  corresponding  to  the  same  initial  data.  Then,  there  is  an 
eq  >  0  and  order-one  constants  Ej,  0  <  j  <  k,  such  that  for  0  <  e  <  Eq  and 
0  <  t  <  e"3/2^ 

K{j)(‘,t)  -  0(j)(*>t)l  5  Ejte^/2  +  9/4,  (61) 

for  0  <  j  <  k,  and,  if  k  >  1, 

sup  IC(j)^x.t)  -  o*j)(x,t)|  <  Ejte^+5^2,  (62) 

x  ®  R 

for  0  <  j  <  k-1. 

COROLLARY  6.  Assume  the  hypotheses  and  notation  of  theorem  1.  Then 
there  Is  an  e0  >  0  and  order-one  constants  Fj,  0  <  j  <  k,  such  that  for 
0  <  e  <  e0  and  0  <  t  <  c"3^2, 

IC(J)(* *t)B  -  FjcJ/2  +  3/4,  (63) 

for  0  <  j  <  k,  and  if  k  >  1, 

sup  lt(j)(‘.t)|  <  FjeJ/2  +  X.  (64) 

x  c  R 

Remark .  Corollaries  5  and  6,  which  are  concerned  with  the  model  equation 
(13),  are  immediate  consequences  of  theorems  2  and  3,  and  their  corollaries, 
once  theorem  1  is  established.  For  example,  to  obtain  (63),  proceed  as 
follows.  For  e  <  eQ,  we  have  from  (41)  and  (27), 


J 
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||C(j)(*.t)||  <  ||C(j)(*.t)  -  n(  j  j  ( *  »t  )||  +  ||  j  ^  (•  »t )  || 

<  m^'2  +  7/4(e3/2t)  ♦  D^'2  +  3/4, 

at  least  for  0  <  t  <  e"3/2.  Take  Ej  *  Dj  Mj,  which  is  an  order-one 
quantity  since  Dj  and  Mj  both  have  this  property.  The  result  then  follows. 

Corollaries  5  and  6  can  be  improved  somewhat  in  the  sense  that  weaker 
hypotheses  concerning  g  suffice  for  the  stated  conclusions.  The  conclusions 
themselves  cannot  be  strengthened  significantly,  as  is  shown  in  the  next 
section.  However,  all  the  results  of  these  corollaries,  as  well  as  the  main 
theorems,  remain  valid  for  temporal  intervals  of  the  form  0  <  t  <  Te"3^2, 
where  T  >  1  is  fixed,  but  arbitrary,  though  the  resulting  order-one  constants 
appearing  in  the  stated  conclusions  are  then  dependent  on  T. 

An  interesting  case  in  corollary  4  is  j  *  0,  when  (60)  yields  the 
estimate 

sup  he(x,t)  -  ce(x,t)|  <  Nne7^2t,  (65) 

holding  for  0  <  t  <  e"3^2.  Inequaltiy  (65)  Is  exactly  that  obtained  In  (9) 
when  comparing  solutions  of  the  two  linearized  model  equations.  Note  In 
particular  that  at  t  *  t^  *  e"3^2, 

sup  lne(xvt1)  -  ce(x,t1)|  <  Nqc2, 

showing  explicitly  that  the  two  solutions  are  the  same  to  the  formal  order  of 
accuracy  e2  achieved  by  either  model.  Recall  from  section  two  that  tj  Is  the 
temporal  scale  over  which  it  Is  expected  that  significant  modifications  of  the 
wave  profile  will  occur,  due  to  the  accumulation  of  nonlinear  and  dispersive 
effects. 


Whilst  tj  is  indeed  a  long  time  scale  for  the  present  considerations,  it 
is  nevertheless  expected  that  the  estimates  (27)  and  (60)  will  continue  to 
hold  for  t  in  the  range  [e”^2,e’^2].  So  far  this  kind  of  result  has  proved 
elusive  to  analytical  methods  for  the  nonlinear  case.  Such  results  are  easily 
established  if  the  nonlinear  terms  In  the  two  model  equations  are  neglected, 
as  in  the  sample  calculation  in  section  two.  Numerically  obtained  evidence 
supports  the  validity  of  the  estimates  (27)  and  (60)  on  the  longer  time  scale 
t2  *  c"5^2.  This  and  some  other  results  are  presented  in  the  next  section. 


5.  F  JRTHER  COMPARISONS  OF  THE  MODEL  EQUATIONS 


Additional  amplification  and  interpretation  of  the  theory  developed  in 
section  four  is  provi  ded  here.  This  is  accomplished  principally  by  way  of 
some  specific  example*-.  In  the  present  discussion,  we  retain  the  notation  and 
scaling  appropriate  to  (26)  in  section  four.  In  particular,  e  continues  to  be 
used  as  the  measure  c*  the  small  a_,olitude  of  the  initial  wave  profiles. 
According  to  the  presentation  in  section  two,  e"1/2  is  therefore  assumed  to  be 
a  measure  of  the  length  scale  characterizing  these  profiles. 

It  is  worth  reit-erating  that  e‘3^  is  the  smallest  time  scale  over  which 
the  nonlinear  and  dispersive  terms  in  either  model  equation,  written  in  the 

forms  (26),  act  effectively  to  alter  the  shape  of  the  initial  wave  profile. 

-3/2 

On  time  scales  significantly  smaller  than  e  ,  theorem  2  and  its  corollary 
show  that  the  dominant  effect  experienced  by  the  initial  profile  with  the 
passage  of  time  is  induced  by  the  top-order  portion,  ot  +  ox  =  0,  of  the 
equations,  and  is  therefore  simply  translational.  As  pointed  out  in  (65), 
solutions  of  the  two  model  equations  corresponding  to  the  same  initial  data, 
as  In  (26),  differ  by  order  at  most  over  the  entire  temporal  Interval 
[0,e“^2].  For  the  (models  written  in  the  form  (26),  both  and  ce  are 
quantities  of  order  e*  Thus  their  difference  Is  seen  to  be  a  factor  of  e 
smaller  than  the  functions  themselves.  This  Is  exactly  the  formal  order  that 
the  terms,  neglected  in  the  derivation  of  these  equations,  would  be  expected 
to  contribute  over  tn's  time  scale.  Put  another  way,  both  ne  and  ce  have 
formal  resolution  of  order  e  over  the  time  interval  [0,e"3/2].  Hence  our 
results  show  unambiguously  that  on  this  time  scale  the  solutions  of  the  two 
equations  are  the  ^ar^to  the  forrvl  order  of  approximation  afforded  by  either 
equation. 


35- 


Two  issues  arise  naturally  upon  further  consideration  of  theorens  1,2, 
and  3.  First,  the  sharpness  of  the  results  deserves  consideration.  Second, 
it  seems  probable  tha.  the  estimates  expressed  in  (27)  and  (60)  are  valid  on 
the  longer  time  scale  t2  =  e~^,  as  asserted  in  our  conjecture  at  the  end  of 
section  two.  For  the  linearized  model  equations  It  may  be  readily  demon¬ 
strated  tnat  (27)  and  (60)  are  sharply  valid  for  t  in  the  range  [0 ,e-^2]. 
Moreover,  as  pointed  out  in  section  two,  formal  arguments  indicate  that  (27) 
and  (60)  are  sharply  valid  for  t  in  [0,e-^]  for  the  nonlinear  problem  as 
well.  Both  of  these  issues  are  addressed  below,  though  neither  has  been 
conclusively  resolved. 

To  fix  ideas  in  the  present  context,  let  us  agree  to  call  an  inequality 
of  the  general  form. 


1  -  Cjef* 


where  Cj  is  an  order-one  constant,  sharp  if  there  is  another  order-one 
constant  C.  such  that 

w 


Here  eE  Is  a  solution,  or  difference  of  two  solutions,  of  equations  (28)  or 
(30),  r  is  some  fixed  real  number  and  i  i  denotes  some  norm,  not  necessarily 
that  of  L2. 

One  way  to  test  the  sharpness  of  the  results  obtained  thus  far  is  by 

resort  to  examples.  A  particularly  simple  class  of  examples  Is  obtained  by 

1/2 

choosing  the  initial  data  eg(e  x)  to  be  a  small -amplitude  solitary-wave 
solution  of  (26a),  namely 

g(x)  *  3sech2(^  x). 


(66) 


This  function  is  an  element  of  h",  so  it  certainly  satisfies  the  hypotheses  in 
tneorems  1,  2,  and  3  for  any  value  of  k  or  m.  For  e  >  0  given,  the  exact 
solution  of  (26a)  for  this  choice  of  g  is 


S£(x,t)  = 


,  1/2 

sech^{^-  [x  - 


(1  +  e )t] ) . 


(67) 


For  this  special  similarity  solution,  it  is  obvious  that  the  norms  appearing 
in  theorems  3  and  corollary  3  do  not  vary  with  t.  Hence,  for  all  j  >  0, 


|ajse(.,t)n  =  iajCeg(E1/2*  )]i 


j/2  +  3/4 


t 


where  dj  is  a  positive  constant  determined  by  g.  Similarly,  for  all  j  >  0, 

sup  I 3^S  (x,t)|  =  sup  l3j[eg(e1/2x)]|  =  f.ej/2  +  *, 
x«  R  *  e  x  «  R  x  J 

where  again  the  fj  are  positive  constants  determined  by  g.  These  simple 

calculations  show  that,  in  general,  the  bounds  obtained  in  (41)  and  (43) 

cannot  be  improved.  We  have  not  proved  a  theorem  of  genericity,  that  all,  or, 

more  likely,  nearly  all  sufficiently  smooth  choices  of  g  result  in  solutions 

that  sharply  bey  (41)  and  (43). 

The  examples  In  (67)  may  also  be  used  to  show  that,  in  general,  theorem  2 
and  corollary  2  are  sharp.  For  the  exact  solution  of  (28)  subject  to  the 
Initial  condition  cgU^x)  is, 

9  1/2 

oe(x,t)  *  3c  sechM— j-  [x  -  t]).  (68) 

i 

The  functions  in  (67)  and  (68)  are  identical  except  that  they  propagate  at 
slightly  different  speeds.  Because  of  this,  they  draw  apart  and 
the  norm  of  their  difference  grows.  A  straightforward  calculation  shows  that 
U4)  and  (38)  are  sharply  verified  in  this  particular  instance.  Thus  both 


theorem  2  and  corollary  2  are  sharp  In  general.  Again,  a  generic  result  along 
these  lines  has  eluded  us. 

Finally,  we  try  our  example  (67)  in  the  context  of  the  initial-value 
problems  for  (A)  and  (B)  posed  in  (26).  There  is  not  available  a  closed- form 
solution  of  (26b)  corresponding  to  the  initial  data  egU^x),  with  g  as 
given  in  (66).  Consequently,  we  have  had  to  rely  upon  a  numerical  integration 
of  (26b).  The  numerical  scheme  used  for  these  experiments  has  been  described 
in  bona,  Pritchard  and  Scott  (1980).  It  results  essentially  from  discretizing 
an  integal  equation  that  is  equivalent  to  (26b).  Ihe  resulting  scheme  is 
quite  efficient,  and  has  been  proved  to  be  unconditionally  stable  and  fourth- 
order  accurate  in  both  the  spatial  and  temporal  mesh  size.  Extensive 
convergence  tests  for  this  scheme  have  been  carried  out  using  the  solitary- 
wave, 

sc(x,t)  »  X  sech2(|  (rLTr)1/Z[x  -  (c  +  l)t]).  (69) 

For  aqy  C  >  0,  this  is  an  exact  solution  of  (26b),  and  the  convergence,  or 
lack  thereof,  of  a  numerical  scheme  may  be  conveniently  tested  on  It. 

In  tables  one  and  two,  we  show  the  errorsarising  tfien  our  scheme  was 
used  to  Integrate  (26b)  (In  the  variables  appurtenant  to  (10b))  with  the 
initial  condition 

h(x)  *  X  sech^Eg^3^  ^]^2x). 

The  errors  plotted  are  with  regard  to  the  discrete  L2  norm  of  the  difference 
between  the  exact  solution 

U(x,t)  -  2C  sech2([7n-3^-n.]1/2[x  -  (C  +  l)t]), 


(70) 
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Ej  =  l  ^uij  "  u( iAX» ^2^x>  1/2» 


and  that  predicted  by  our  scheme.  If  u^j  is  the  value  given  by  the  numerical 
scheme  at  the  point  (iAx,jAt),  where  ax  and  At  are,  respectively,  the  spatial 
and  temporal  mesh  size,  then  the  error  at  the  jth  time  step  is  defined  to  be 

(71) 

i 

and  the  relative  error  at  the  jth  time  step  is  taken  to  be 

ei  =  - ^ - 2 - T77'  (72) 

J  (£  UfiAx.jAtrAx}1^ 

i 

It  is  the  relative  errors  that  are  tabulated.  (Note  that  the  summation  in  the 
above  formulas  must  be  truncated.  Because  of  the  rapid  spatial  decay  of  the 
solution  in  question,  such  a  truncation  can  be  made  without  sensibly  affecting 
the  approximation  to  the  L2  norm.  In  all  cases,  spatially-truncated  values  of 
u^j/2C  were  less  than  10“8  in  our  computations.)  Preliminary  experiments  had 
shown  that  the  choice  Ax  ssrAt  gave  the  best  results  (in  the  sense  of  accuracy 


achieved  for  the  work  expended)  and  for  these  computations  we  have  taken 
Ax  =  At  =  A  ,  say. 

Table  1.  The  relative  error  e«  induced  in  integrating  a  solitary 
wave  (70)  of  amplitude  2C  =  1.0 

(ax  *  At  =  a;  an  entry  in  a  row  labelled  'ratio'  is  the  ratio  of  the 
numbers  above  and  below  that  entry.) 


a  Time 

0.640 

9.600 

30.720 

72.320 

120.320 

c 

0.16 

0.147E-3 

0.688E-3 

0.249E-2 

0.141E-1 

0.400E-1 

ratio 

15.7 

17.2 

28.4 

31.3 

31.3 

r 

0.08 

0.937E-5 

0.400E-4 

0.876E-4 

0.451E-3 

0.128E-2 

ratio 

15.9 

16.7 

26.5 

32.2 

32.4 

0.04 

0.588E-6 

0.240E-5 

0.330E-5 

0.140E-4 

0.395E-4 

ratio 

15.7 

16.1 

21.4 

24.4 

28.8 

0.374E-7 

0.149E-6 

0.154E-6 

0.574E-6 

0.137E-5 

Table  2.  The  relative  error  e*  induced  in  integrating 
a  solitary  wave  (70)  of  amplitude  2C  =  0.1 


(ax  =  At  =  a;  an  entry  in  a  row  labelled  'ratio'  is  the  ratio 
of  the  number  above  and  below  that  entry.) 


A  Time 

0.640 

9.600 

30.720 

72.320 

120.320 

0.32 

0.227E-3 

0.334E-2 

0.108E-1 

0.256E-1 

0.429E-1 

ratio 

16.0 

15.7 

15.8 

15.7 

15.7 

0.16 

0.142E-4 

0.213E-3 

0.685E-3 

0.163E-2 

0.274E-2 

ratio 

16.0 

15.9 

15.9 

15.8 

15.9 

0.08 

0.890E-6 

0.134E-4 

0.430E-4 

0.103E-3 

0.172E-3 

ratio 

16.0 

16.0 

16.0 

16.1 

15.9 

0.04 

0.557E-7 

0.837E-6 

0.269E-5 

0.641E-5 

0.108E-4 

Having  determined  the  sort  of  accuracy  inherent  in  the  numerical  scheme, 
we  then  introduced  the  initial  data  eg(e^2x),  where  g(x)  *  2  $ech2([|ol^2x) 
(corresponding  to  (66),  except  for  the  equation  in  the  form  (10a)).  Equation 
(10b)  was  integrated  numerically,  with  the  just -mentioned  form  of  initial 
data,  for  a  range  of  values  of  e.  At  time  t,  the  difference  between  the 
numerically  computed  solution  of  (10b)  and  the  exact  solution  of  (10a)  was 
formed.  Let  M£(t)  denote  the  maximum  value  of  this  difference,  and  consider 
the  function  log(Me(t)).  Because  of  the  estimate  (65),  it  is  expected  that 

1 og(Me (t ) )  2  constant  +  g  log  e  +  log  t. 

Attention  is  fixed  on  the  particular  times,  tQ  *  e”^2,  tj  *  e“3^2,  t2  s  e~5 
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and  t3  *  e  considered  earlier.  For  these  times,  we  expect, 

log(M£(tj))  2  constant  +  (7  -  7  -  j >1 og  e  -  constant  +  (3  -  j)log  e. 


An  idea  of  how  well  this  relation  is  obeyed  may  therefore  be  obtained  by 


plotting  log  He(tj)  versus  log  e  for  various  small  values  of  e .  This  is 
shown  in  figure  1.  The  general  pattern  of  the  results  appears  to  confirm  that 
(65)  is  sharply  valid  over  the  entire  temporal  range  [Q,e"^3,  at  least  for 
this  particular  example.  In  table  3,  the  results  plotted  in  figure  1  are 
tabulated,  along  with  the  numerical  values  of  the  slopes  determined  therefrom. 


Remark.  A  direct  comparison  of  (A)  and  (B)  could  also  be  made  by  posing 
as  initial  data  a  small -amplitude  solitary -wave  solution  of  (B).  One  could 
then  attempt  to  use  the  inverse-scattering  theory  to  infer  properties  of  the 
resulting  solution  of  the  KdV  equation  (A).  We  have  investigated  this  possi¬ 
bility  and  have  found  the  fruits  of  our  labour  to  be  suggestive,  but  not 
conclusive.  For  the  expert,  it  is  worth  adding  that  we  were  not  able  to  gain 
sufficient  control  of  the  dispersive  tail  to  effect  the  needed  estimates. 

A  final  and  more  delicate  way  in  which  these  two  model  equations  were 
compared  Is  now  described. 

It  is  an  interesting  fact  that  certain  classes  of  solutions  of  the  KdV 
equation  (A)  have  a  rather  simple  structure  for  large  values  of  t.  Basically, 
the  asymptotic  form  of  these  solutions  comprises  a  widely-spaced  sequence  of 
independently  propagating  solitary  waves  followed  by  a  small -amplitude 
dispersive  tail.  This  behaviour  was  first  noted  numerically,  and  was  then 
established  analytically  by  use  of  the  so-called  inverse-scattering  theory  for 
KdV  (see,  for  example,  Mlura  1974,  1976). 

It  is  natural  to  inquire  whether  or  not  the  same  property  holds  for 
equation  (B).  A  related  property  of  the  KdV  equation,  whereby  solitary  waves 
recover  their  exact  form  upon  emerging  from  an  interaction,  appears  to  be 
false  for  (B)  (see  Bona  1980  and  the  references  included  therein). 

This  latter  result  does  not  preclude  the  possibility  that  solutions  of  (B) 
exhibit  the  aforementioned  asymptotic  form.  The  issue  seems  to  be  beyond  the 
reach  of  the  analytical  tools  currently  at  our  disposal.  Consequently,  it  has 
been  investigated  using  the  numerical  Integration  procedure  for  (B)  discussed 
earlier. 


We  present  here  an  example  of  the  sort  of  outcome  observed  in  a  number  of 
numerical  experiments.  In  this  example,  the  initial  data  was  taken  to  be 

c(x,0)  =  exp(-x2/10).  (73) 

The  calculations  were  performed  with  ax  s  at  *  0.16,  and  a  check  was  made  by 
comparing  with  the  same  integration  where  ax  *  At  *  0.08.  The  procedure 
appeared  to  have  converged,  for  all  practical  considerations.  As  always,  the 
calculations  were  run  on  a  finite  spatial  grid  which  was  chosen  to  be  large 
enough  that  values  of  c  truncated  by  this  limitation  did  not  exceed  10“®.  A 
feature  of  our  numerical  procedure  was  a  peak-finding  subroutine  which,  at 
each  discrete  time  step,  located  local  maxima  in  the  bulk  of  the  solution  and 
compared  the  wave  profile  near  such  a  peak  with  a  solitary -wave  solution  of 
(10b)  having  the  corresponding  amplitude  and  positioned  so  that  its  peak 
coincided  with  the  local  maximum  In  question.  (This  routine  ignored  peaks  of 
amplitude  less  than  0.01.)  The  routine  also  computed  the  local  speed  of  the 
wave  based  on  the  movement  of  Its  crest.  The  result  of  our  calculation  is 
pictured  In  figure  2,  which  Is  now  described  In  some  detail.  All  differences 
between  the  computed  profile  and  various  solitary-wave  solutions  of  (10b)  are 
relative  discrete  L2-norms,  as  defined  In  (72). 

1)  The  Initial  profile  differed  from  a  solitary  wave  of  the  same 
amplitude  by  0.963. 

2)  Two  peaks  were  visible  by  time  6.4,  one  of  amplitude  1.307  and  another 
of  amplitude  0.4442. 

3)  By  time  16.0,  three  peaks  had  resolved  themselves,  with  respective 
amplitudes  1.358,  0.4633  and  0.0883. 

4)  By  time  35.2,  four  peaks  were  discernible  with  amplitudes  1.357, 
0.4738,  0.07206,  and  0.01072. 

5)  A  fifth  peak  had  emerged  by  time  105.6.  The  amplitudes  of  the  peaks 
at  this  time  were  1.3&1,  0.4748,  0.0595,  0.0161,  0.0101. 


7 5.04 


The  development  over  time  of  the  numerical .approximation 
solution  of  ( 10b)  with  Inltal  data  exp(-x‘/10). 
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Now  the  finer  structure  of  the  evolution  of  the  various  peaks  is  des¬ 
cribed.  In  the  ensuing  discussion  the  terms  'speed  error'  and  'shape  error' 
are  used  respectively  to  mean  the  difference  in  speed  and  shape  between  the 
computed  waveform  and  a  solitary  wave  (70)  having  the  same  amplitude  2C  as 
that  of  the  computed  wave. 

Peak  no.l.  The  first  peak  grew  steadily  in  amplitude  to  a  maximum  value 
of  1.358  at  about  t  =  13.  Thereafter  its  amplitude  decreased  very  slowly  to  a 
value  of  1.348  at  t  =  144.  This  latter  decrease  is  believed  to  result  from 
numerical  errors;  such  slow  attrition  was  consistent  with  our  Integration  of 
solitary-wave  solutions  of  (10b). 

At  t  *  12.8,  when  the  wave  had  virtually  reached  its  ultimate  height,  its 
speed  of  propagation  was  1.679,  which  differed  by  0.3E-2  from  that  of  a 
solitary -wave  solution  of  (10b)  with  the  same  amplitude.  The  difference 
between  the  profile  near  this  first  crest  and  a  solitary  wave  of  the  same 
height  was  .067.  At  t  *  32.0,  the  speed  differed  from  that  of  the  "appropriate" 
solitary  wave  by  only  0.35E-3  and  the  error  In  shape  was  0.89E-3.  At  t  =  48.0 
and  144.0  the  speed  error  was  .35E-3  and  .34E-3,  respectively,  and  the  shape 
error  was  .89E-3  and  .89E-3,  respectively. 

Peak  no.  2.  This  crest  emerged  fairly  quickly,  beginning  its  independent 
existence  with  an  amplitude  of  0.4480,  and  growing  steadily  in  amplitude  to  a 
maximum  of  0.4748  at  about  time  55.  It  held  this  latter  value  thenceforth. 

At  t  *  9.6,  the  speed  of  this  wave  as  determined  by  the  movement  o^  the  peak 
was  1.208,  which  is  0.015  slower  than  a  solitary  wave  of  Its  amplitude.  The 
difference  In  shape  was  2.476. 

Both  the  error  In  speed  and  shape  decreased  rapidly  as  this  peak  became 
more  Isolated.  At  t  =  16.0  and  32.0  the  speed  errors  were  .36E-2  and  .96E-3, 
respectively,  and  the  shape  errors  were  .72  and  .072,  respectively. 
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When  the  amplitude  had  stabilized,  at  about  t  =  55,  the  speed  error  was 
.42E-3  whilst  the  shape  error  was  .68E-2.  At  t  =  144.0,  the  speed  difference 
was  .37E-3  and  the  shape  difference  .14E-2. 

Peak  no.  3.  This  peak  first  emerged  at  time  about  16  with  an  amplitude 
of  0.08828.  Its  speed  at  this  time  was  0.9477  and  it  differed  markedly  from  a 
solitary  wave  (shape  error  5.73).  The  amplitude  decreased  monotonically  taking 
the  values  0.07339,  0.06458,  0.06038  and  0.05695  at  t  =  32.0,  64.0,  96.0  and 
144.0,  respectively. 

The  speed  of  the  crest  increased  monotonically,  surpassing  1  by  t  *  28.8. 
At  this  time  its  speed  differed  from  that  of  a  solitary  wave  by  0.0345  and  the 
shape  difference  was  4.13.  These  two  measures  continued  their  development  as 
f ol 1 ows . 

At  t  =  48.0  the  speed  error  was  0.165E-1  and  shape  error  was  4.51 
64.0  "  0.117E-1  "  4.66 

96.0  "  0.708E-2  "  0.50 

128.0  "  0.483E-2  "  0.46 

144.0  "  0.417E-2  “  0.44 

Peak  no.  4..  This  peak  first  emerged  at  t  «  35  with  an  amplitude  of 
0.01072.  Its  initial  speed  and  shape  were  quite  different  from  those  of  the 
appropriate  solitary  wave;  the  relevant  errors  were  0.1129  and  43.8, 
respectively.  The  amplitude  of  the  crest  initially  increased  to  a  value  of 
0.01421  at  t  !  51.2,  to  0.01540  at  t  =  64.0  and  reached  a  maximum  of  0.01615 
at  t  *  96.0.  The  amplitude  thereupon  decreased  slowly,  taking  values  0.01586 
at  t  =  128.0  and  0.01560  at  t  =  144.0. 

The  speed  of  this  crest  Increased  with  time,  taking  a  value  of  0.9091  at 
t  *  51.2,  a  value  of  0.9221  at  t  *  64.0,  and  a  value  of  0.9417  at  t  *  96.0. 


For  t  in  excess  of  96.0  the  speed  of  the  crest  continued  to  increase, 
even  though  the  amplitude  was  now  decreasing,  taking  the  values  0.9527  and 
0.9564  at  t  =  128.0  and  144.0,  respectively.  At  t  =  144.0,  this  wavelet 
differed  in  shape  from  that  of  a  solitary  wave  by  2.96. 

Peak  no.  5.  This  peak  first  emerged  near  t  =  105.6  with  an  amplitude  of 
0.01008.  It  seemed  to  have  grown  from  zero.  Its  amplitude  then  increased 
with  time  to  a  value  of  0.01097  at  t  =  144.0. 

The  initial  speed  of  the  crest  was  0.9135  at  t  =  108.8.  This  rose 
steadily  to  the  value  0.9291  at  t  =  144.0. 

As  the  reader  will  easily  discern  from  the  foregoing  description,  it 
would  be  optimistic  to  claim  that  the  full  picture  of  the  evolution  of  the 
Initial  data  (73)  under  the  action  of  (10b)  is  captured  by  our  numerical 
calculations.  Fully  realizing  the  need  for  some  caution  here,  we  nevertheless 
feel  the  following  summary  Is  an  accurate  description  of  what  really  happens 
to  this  Initial  data. 

It  appears  that  three  solitary  waves  have  emerged  from  this  Initial  pro¬ 
file,  though  the  third  solitary  wave  still  had  considerable  evolution  to 
undergo  before  It  could  reasonably  be  said  to  have  established  Its  asymptotic 
form.  The  fourth  and  fifth  waves  were  surely  part  of  a  ‘dispersive  tail'.  They 
showed  no  signs  of  settling  down  to  a  uniform  amplitude.  More  Importantly 
their  speeds  were  significantly  less  than  1.  By  consulting  the  formula  (70), 
one  determines  that  a  solitary-wave  solution  of  (10b)  propagates  with  a  speed 
exceeding  1.  Moreover,  the  other  details  of  these  last  two  waves  fit  within 
the  general  structure  observed  for  dispersive  tails  arising  In  the  integration 
of  (10b)  (see  Bona  et_al_.  1980). 

If  one  takes  the  view  that  this,  and  other  like  calculations,  do  point  to 
(10b)  having  the  property  that  certain  initial  profiles  resolve  themselves 


into  a  sequence  of  solitary  waves  and  a  dispersive  tail,  then  a  host  of 
questions  present  themselves.  These  questions  will  not  be  addressed  here. 

But  it  is  worth  pointing  out  that  if  indeed  our  surmise  is  correct,  then 
solutions  of  the  two  model  equations  corresponding  to  physically  relevant 
initial  data  may  agree  qualitatively  over  indefinitely  large  time  scales. 

This  would  be  the  case,  for  example,  if  by  the  time  t2  =  e"^2,  both  solutions 
had  already  sorted  themselves  into  more  or  less  independently  propagating 
solitary  waves.  For  then,  by  theorem  1,  these  solitary  waves  rust  correspond 
to  each  other  in  number  and  be  very  close  In  amplitude.  The  further  evolution 
of  these  solutions  will  then  be  quite  simple;  an  observer  far  downstream  would 
report  of  both,  seeing  like  sequences  of  solitary  waves  trailed  by  a  small - 
amplitude  dispersing  disturbance. 
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